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Abstract. It is "well known" that there is no explicit expression for the 
Black-Scholes implied volatility. We prove that, as a function of underlying, 
strike, and call price, implied volatility does not belong to the class of in- 
finite functions. This does not rule out all explicit expressions, but shows that 
implied volatility does not belong to a certain large class, which contains many 
elementary functions and classical special functions. 



1. Introduction 

The Black-Scholes formula is the most popular way to price European stock 
options. Assuming a frictionless market and an underlying modelled by geometric 
Brownian motion with volatility a, the arbitrage free price of a European call option 
with strike K, maturity T, and initial stock price S is 

C BS (S,K,T,a) = SN(dx) - KN(d 2 ), 
where iV" is the cdf of the standard normal distribution, and 



\og(S/K) aVT 

"1,2 = 7= ± 



<jVT ' 2 

(We assume a zero interest rate for simplicity.) Given a call price surface C that 
satisfies the no-arbitrage bounds 

(1) (S- K)+ < C(S,K,T) < S, 

the Black-Scholes formula can be inverted w.r.t. a to obtain the implied volatility 

f^imp ■ 

(2) C BS (S,K,T,a iiap (S,K,T)) = C(S,K,T). 

Implied volatility allows to compare option prices for different strikes, maturities, 
underlyings, and valuation times. For entry points to the literature on implied 
volatility, see, e.g., EH HU HI HP • 

The compute the implied volatility from a call price, equation © has to be 
solved numerically; many authors have added that the inversion cannot be done 
in closed form. The latter statement has been made so often that one might ask 
whether it can be turned into something tangible. First, note that it is hard to say 
anything about explicit expressions for the function 

(3) (S,K,T)^a imp (S,K,T). 

This function actually depends on C(-,-,-), and it would be unnatural to make 
assumptions on the generic call price surface C that would allow to refute a closed 
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form for (J3J. Instead, we consider the function I which satisfies 

(4) C BS (S, K, 7(5, K, c)) = c, S, K > 0, (S - X)+ < c < S, 
and is defined on the open set 

Z>j := {(S, K, c) e M 3 : 5, If > 0, (S - K)+ < c < S}. 

The maturity T has been omitted here, and we will continue to do so, assuming 
a fixed T > throughout. Note that any kind of explicit expression for 7 with 
variable T would imply an explicit expression for any fixed T . The function I 
clearly exists and is unique, since the Black-Scholes call price increases from the 
lower to the upper no-arbitrage bound as a increases. If a call price surface C 
satisfies the slightly strengthened bounds 

(5) (S-K)+ < C(S,K) < S, 

which usually hold for T > 0, then the corresponding implied volatility <7i mp satisfies 

a imp (S,K) = I(S,K,C(S,K)), S,K>0. 

An example where the lower bound in ((SJ) fails would be an out-of-the-money call 
that cannot move into-the-money, as may happen in a binomial model, say. We 
impose (J5J instead of (JTJ to make the domain of I open, so that I is differentiable 
on its whole domain T>i. 

Since the Black-Scholes call price is real analytic for S,K,o~ > 0, and the Black- 
Scholes vega dC-Qs/da is positive, the implicit function theorem [15] shows that I is 
real analytic (hence C°°) on T>j. The question we now ask is whether the function I 
admits a closed form. To give a partial answer, we recall in Section [2] the definition 
of the class of D-finite functions. The main result of this note (Theorem [1] in 
Section [3] below) is that I is not D- finite. 

We briefly mention some other results of the kind "a certain function does not 
belong to a certain class". There is a wealth of theorems about transcendental 
functions, i.e., non-algebraic ones; for examples, see, e.g., Schmidt [3D] and Bell et 
al. [3] and the references given there. Rubel [TSJ discusses "transcendentally tran- 
scendental" functions, which means that they do not satisfy an algebraic differential 
equation. The class of D-finite functions, which we are interested in, lies in between 
the algebraic and differentially algebraic classes. Finally, Bronstein et al. [S] have 
shown that the Lambert W function is not Liouvillian, which roughly means that 
it cannot be expressed by iterated integration and exponentiation starting with an 
algebraic function. 

2. D-FINITE FUNCTIONS 

Suppose that a C°°-smooth function / is defined on an open set T>f C W 1 . It is 
called D-finite if it satisfies PDEs 

+ iU-i^ TFT + • • ■ + ^uWf^ + ^i,o(x)/(x) = 0, 
ox^ axi ox\ 

(6) ; 

P„, dn (x)^-^ + P w , dn _ 1 (x) C; +■■■ + P na ( x )^l + P„, (x)/(x) = 0, 

ox n n ox n n ox n 

valid for x = (xi, . . . ,x n ) £ T>f, where di > 1 for i = 1, . . . , n, and the Py are 
polynomials such that P,^ is not identically zero on £>/ for i = 1, . . . , n. If / 
is real analytic, let us fix an arbitrary point Xo € T>f and consider the Taylor 
expansion of / at Xo- When viewed as a formal power series, the PDEs © show 
that its partial derivatives generate a finite dimensional vector space over the field 
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of rational functions. This algebraic definition of D-finiteness is the usual one in 
the literature jTTJ [22] . It cannot be used directly for functions, though: Unlike the 
held of formal Laurent series, C°°(T>f) is not a vector space over the field of rational 
functions. Recall that in one dimension, the field of formal Laurent series R((X)) 
in the indeterminate X consists of the formal series 

| a n X n : hq £ Z, a n € K, n > tiq|. 

n>no 

The multivariate case is more involved; see Aparicio Monforte and Kauers [2J. 

Rational, and, more generally, algebraic functions are D-finite, as are the elemen- 
tary functions exp, log, sin, and cos (but not tan and cot) . About 60% of the spe- 
cial functions covered in the classical handbook by Abramowitz and Stegun [I] are 
D-finite. Examples include Bessel functions, Airy functions, exponential integral, 
dilogarithm, and hypergeometric series. Furthermore, the class of D-finite func- 
tions is closed under addition, multiplication, (in-)definite integration, and Laplace 
transform. Division, on the other hand, does not preserve D-finiteness in general, 
nor does composition, unless the inner function is algebraic. Thus, x — > e^ is 
D-finite, whereas x — ¥ e e and x — > l/sina; are not. The closure under algebraic 
substitution will be applied in the proof of our Theorem [T] below. 

D-finite functions resp. power series have been studied extensively in the discrete 
mathematics and symbolic computation literature. Indeed, provided the coefficients 
of the polynomials in ([6]) and the initial conditions can be represented with a 
finite amount of information, one has a finite data structure that can represent 
reasonably general functions (or, in symbolic computation, rather formal power 
series). The closure properties mentioned above are effective, in the sense that 
there is an algorithm that computes the D-finite specification for a sum of two given 
D-finite functions, etc. There are also algorithms for proving identities involving in- 
finite functions, which are about to render tables of such formulas (certain definite 
integral evaluations, e.g.) obsolete to some extent. For further information, see, 
e.g., Chyzak and Salvy [B] or Koutschan |14) . 

Summarizing, the class of D-finite functions is extensive, but still there are plenty 
of "explicit" functions that are not D-finite. The following theorem thus gives a 
partial answer to the question asked in the title of the present note. 

3. Main result 

Theorem 1. The function I : T>j Cl 3 -) (0, oo) defined by ^ is not D-finite. 

Proof. It is often hard to show directly that a multivariate function is not D-finite, 
whereas there are several methods available that deal with univariate functions [HO 
112] . Thus, our strategy is to find a useful specialization of / to a univariate function. 
We let the underlying's initial price be proportional to the strike: S — eK. For 
< K < I/e, the inequality 

(e - l)K < c(K) := (e - l)K + eK 2 < eK 

holds, and so the function 

}{K) :=I(eK,K,c(K)) 

is well-defined for < K < I/e. We assume that I is D-finite, and want to show that 
then / would be D-finite, too, by studying its Taylor series at K = l/(2e). First, wc 
expand the analytic function I in a neighborhood of the point {\,j^,\ — ^)^T^i- 

(7) I(S,K,c)= £ lijkiS -ly(K-±y(c-(i-i;)) k . 

i,j,k>0 
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Since / is L>-finite, so is the formal power series 

Y, lvkX l YiZ k 

i,j,k>0 

in the indeterminates X, Y, Z. (Recall that this means that its partial derivatives 
generate a finite dimensional subspace of the space of formal Laurent series, which 
is a vector space over the field of rational functions M.(X, Y, Z).) The algebraic 
substitution X -> el, Y -t X, Z — s> (e - l)X + eX 2 preserves D-finiteness QH [T71 
122] . and so the univariate formal power series 

(8) Yl (eXyX k (eX + eX 2 ) k 

ij,k>0 

is Infinite. This series represents the analytic function x H i f[x + ^) in a neigh- 
borhood of zero. Indeed, by ([7]), 

i,j,k>0 

and this is (0 with X replaced by K — l/(2e). This shows that £>-finiteness of / 
implies D-finiteness of /. 

We will employ the following variant of the Black-Scholes formula (see Roper 
and Rutkowski [TS]; note that the maturity T > is fixed throughout): 

(9) C BS (S,K,a) = (S - AT + s£^N>(^^l + V -)dv. 
Define the function F : (0,oo) — > (0, 1/e) by 

(10) F(x) :— f ' N '(- + V -)dv. 



v 2, 

Then, for S = eK, the log-moneyness is log(S'/A") = 1, and equation ^ becomes 

C BS {eK,K,a) = (e - l)K + eKF(crVT). 
By the definition of /, we infer that the equation 

(e - l)K + eKF(VTf(K)) = c{K), < K < 1/e, 
must hold, and hence 

F{ VT fm) = = K, 

or 

f(K) = F-\K)/VT, 0<K<l/e. 

The following lemma shows that F~ l is not D-finite, hence / and / are not D-finite, 
either. □ 

Lemma 2. Let the function F : (0,oo) — > (0, 1/e) be defined by (flU)) . Then the 
inverse function F^ 1 : (0, 1/e) — > (0,oo) is not D-finite. 

Proof. We will show that F^ 1 has a singularity at zero whose type is incompatible 
with D-finiteness. First, we determine the asymptotics of F itself as x — > + . For 
< v < x, we have 

l> e -" 2 / 8 > e -* 2 / 8 = l + 0(x 2 ), 

and so 

(11) N'(- + -) = ^=e- 1 ^ 2v2 \l + 0(x 2 )), x^0+. 

2/ \/2e7r 
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Substitution, integration by parts, and the expansion of N at infinity yield 

1 r x i r 00 — v 1 12 

1 f e -l/W dv= J=[ tlL d y 



2?r Jo V2tt J 1/x y 



2tt V2tt 

x 



1/x 

er 1 /^ -l + N(l/x) 

(12) = _^ e -i/(^) ( i + o (a; 2 )). 
From (fTTI) and (fT2)) we conclude 

3 

(13) F(i) = -^=e~ 1/{2x2) {l + 0{x 2 )), x^0+, 

V2e7r 

and thus, taking logarithms, 

(14) _i ogj p( a ,) = _L + (logi), x^0+. 

Now we study the asymptotics of i* 1-1 . For small a; > 0, formula ([13]) yields 
F(x) < e" 1 ^ 2 ^ ). Replacing x by F _1 (x) in this inequality and rearranging, we 
obtain 

(15) F-\x) < - 1 

V21og(l/x) 

Now we replace a; by i 7 '~ 1 (a;) in ([TJ]) and use (|T5"j) in the O-estimate: 

(16) l0 ^-^W + °( l0 S 10 ^ 



This implies that ([15]) is sharp: 

(17) F _1 (.t) 1 x^0+. 

V21og(l/x)' 

Observe that a univariate infinite function g solves an ODE 

P d {x)g {d \x) + ■■■ + Pi(x)g'(x) + P (x)g(x) = 

with polynomial coefficients, where Pa is not identically zero. Such a function has 
an analytic continuation throughout the complex plane, up to a finite number of 
poles and branch points. At each such singularity xq, an asymptotic expansion 
holds, which may feature powers of x — xq, exponential terms, and integral powers 
of log(a; — Xo). This 19th century result |13) is a powerful method for proving non- 
£>-finiteness; see Flajolet et al. [7J [5] for details. In particular, a fractional power of 
a logarithm, as in (fT7j) . is a "forbidden" asymptotic element, and shows that F^ 1 
is not D-fmite. □ 
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